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Abstract 



We explore the relationship between two figures of merit for an adi- 
abatic quantum computation process: the success probability and the 
minimum gap between the ground and first excited states. We study a 
generic adiabatic algorithm and show that a rich structure exists in the 
distribution of these two important variables. In the case of two qubits, 
F^ the success probability is to a good approximation a function of the min- 

^H imum gap, the stage in the evolution at which the minimum occurs and 

>— ' the computation time. This structure persists in examples of larger sys- 

tems and appears to be closely linked to the choice of initial and problem 
K^ Hamiltonians. 

^ 1 Introduction 

1^ The promise of a qualitative advantage of quantum computers over classical ones 

^-v in solving certain classes of problems has led to a massive effort in theoretical 

^^ and experimental investigation of controlled, quantum-coherent systems. The 

^-H standard model of quantum computing is analogous to classical computing in 

^ that it involves the precise manipulation of individual qubits in a register to ap- 

• ^ ply logic gate operations. However, the requirement of precise time-dependent 

^^ control of individual qubits has been found to be hard to achieve experimen- 

^ tally while still maintaining the quantum coherence of the system. A number 

of alternative approaches have been proposed, of which adiabatic quantum com- 
puting (AQC) is a promising example. It involves the evolution of a quantum 
system from a simple Hamiltonian with an easily-prepared ground state to a 
Hamiltonian that encodes the problem to be solved. If the system is initially 
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prepared in the ground state and the time evolution occurs slowly enough to 
satisfy the adiabatic theorem, the final ground state that satisfies the problem 
at hand can be read out with a high probability [1 . 

The AQC process requires a Hamiltonian that interpolates from a simple 
initial configuration, Hi, to one that encodes the problem under consideration, 
Hf. A Hamiltonian for linear interpolation can be written 

H{s) = {l-s)Hi^sHf, (1) 

where s G [0, 1] is the reduced time s = t/T, T being the computation time. 
The instantaneous eigenvalues and eigenstates of the Hamiltonian of an n-qubit 
system are given by 

H{s) \£; s) = Ee{s) \£; s) , with Eo{s) ^ Ei{s) ^ • • • ^ ^2^-1(5). (2) 

The instantaneous state of the system is given by |'0(s)), the solution of Schrodinger's 
equation, which in the reduced time reads 

^^\i,{s)) = -iTH{s)ms)). (3) 

The system is initially in the ground state of Hi'. |?/;(0)) = |0; 5 = 0). 

At the end of the evolution we require a measure of how closely the state 
vector, |'0(1)), corresponds to the desired result, |0;s = 1). This is provided by 
the success probability 

P=\{0;s=imi))\\ (4) 

It was shown in [l] that for P to be arbitrarily close to 1, the following condition 
must be satisfied: 

^»A^ (5) 

min 



where 

t = max ( 1 ; 5 

0<s<l 



0;s 



(6) 



ds 
is of the order of a typical eigenvalue of H and 

{E^{s) - Eo{s)) (7) 
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is the minimum gap, which occurs at reduced time(s) 5* : £^1(5*) — Eo{s'') = 
Amin- If £ is considered constant, by ([5| A^in is the variable which determines 
the T that is required. It is clear that P will also be directly related to T and 
Amin- These two variables, P and Amin, are used interchangeably throughout 
the literature to quantify the performance of a given computation (e.g. contrast 
[2] with [3]), and are assumed to increase monotonically with each other. The 
question of how either of these variables varies with system size, n, is an impor- 
tant one that is often addressed. However, the exact nature of the relationship 
between these two important figures of merit has not been adequately explored. 
We explore the relationship between P and Amin by looking at the statistical 
distributions of these two variables over a large set of generic problem Hamilto- 
nians (Hf). We start by considering a simple two-qubit system and show that 
a rich structure arises in the scatter plots of success probability against Amin- 
We then go on to look at the scatter plots in three-, four- and five-qubit systems 
and find that, although some of the finer details of the structure are washed out, 
some remain. We propose that the structure in the distributions arises through 
our choice of Hamiltonians. 



2 A generic adiabatic algorithm 

We wish to look at the distribution of the success probabihty and Amin over a 
large set of problem Hamiltonians. We use a simple, yet generic, model that is 
scalable and can be readily solved numerically. For iJ^, we use a Hamiltonian 
that is easy to construct and whose ground state can be readily found: 

n n 

i:^i = - ^ cr^'^ = - ^ I (g) • • • (g) cr^ (g) • • • (g) I, (8) 

i=l i=l 

where cFa^ with a = x^ y oi z^ denote the usual Pauli matrices, n is the number 
of qubits in the system and the index i denotes which qubit the operator is 
applied to. This Hi is simply a transverse field acting on all the qubits and its 
ground state is an equal superposition of all 2^ computational basis states. 

For Hf^ we use a Hamiltonian where all possible couplings between the n 
qubits are realised in the z-direction, with random strengths: 

where Xi is the ith digit in the binary representation of x and Jx are the cou- 
pling coefficients. These will be selected from a given random distribution. Hf 
is diagonal in the computational basis so that the binary-ordered set of states 
\y) is a permutation of the energy-ordered set of states \l;s = 1) defined in eq 
[2J A Hamiltonian of this type can easily be used to encode any finite compu- 
tational optimisation problem (minimisation of a function / : {0, 1}^ ^ M) by 
choice of the {Jx}^ including NP-hard problems such as the travelling salesman 
problem (TSP). It is important to note that, at present, only one- and two-qubit 
interactions are experimentally feasible, although the TSP can be implemented 
with this restriction [4 . 

For each sample in the scatter plots, we solve the Schrodinger equation 
numerically over the parameterised time, 5, for a given computation time, T, 
using the Dormand-Prince method [5 . This is an adaptive step-size algorithm; 
solutions accurate to fourth- and fifth-order in As are used to estimate the local 
error in the former. If it is less than the desired tolerance, then the fifth-order 
solution is used for the integration. Otherwise the step-size As is decreased. 

3 Two-qubit simulations 

Fig. [l] is a scatter plot of success probability against minimum gap for a large 
set of problem instances, with the coupling coefficients drawn from the uniform 
distribution ^/(— 3, 3). The computation time is T = 5. Observe the sharp upper 
and lower edges. The lower bound of the success probability is always 1/4 for 
inffnitesimally small Amin- This arises when Ji = J2 = J3 = 0, which means 
there is four-fold degeneracy at 5 = 1 and the system remains in its original 
ground state. 

It is important to verify that this structure is independent of our choice of 
random distribution of coupling coefficients and that it is also not an artefact 
of the pseudo-random number generators used. Fig. [2] also shows scatter plots 
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Figure 1: Success probability against minimum gap for a two-qubit system at 
a computation time of T = 5. The Ji have been chosen from the uniform 
distribution Z//(— 3,3) for the 500,000 random problem instances. The data 
points are coloured by | J3I. 



of success probability and minimum gap, but in this case the coupling coeffi- 
cients are drawn from a Gaussian distribution, A/'(0, 1^) (mean /i = 0, standard 
deviation <j = 1). The trends and structure in the distributions are similar to 
those shown in Fig. [l] However, there are some subtle differences in sharpness 
between the Gaussian and uniform cases. For a large minimum gap, the lowest 
probability occurs for large IJ3I, so we see a sharp cutoff in the uniform case 
and a rougher, more sparsely populated edge in the Gaussian case. In general 
though, this shows that the results are independent of our choice of coupling 
coefficients and, as a different psuedo-random number generator routine was 
used, we can say that the results are not a numerical artefact. 

Four computation times are shown: T = 5, 10, 20 and 40. As this increases, 
the distribution shifts and tends towards a success probability of 1, as expected 
from the adiabatic theorem. 

The two interesting features of these scatter plots are the well-defined sharp 
edges and the densely-populated bands. It is clear that the bands correspond 
to groups of Hamiltonians with similar | J3I. The bands where J3 = can be 
seen as two separable one-qubit evolutions for Ji and J2, so the total success 
probability is simply the product of the one-qubit success probabilities. Another 
interesting point to note is that the bands of similar Hf gradually reverse in 
order in the distribution as the computation time T is changed. 

We have supplemented the uniform random data with sets of Ji chosen on 
a rectangular grid with the same cut-offs. These have the advantage that all 
problem Hamiltonians with a given value of J3 can be plotted in the (Ji, J2)- 
plane and coloured by their minimum gap or success probability; see Fig. |3] 
These two plots can be considered ffist- and second-order stability diagrams 
respectively. 

The ffist part of Fig. |3] (minimum gap) can be explained as follows: dark 
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(a) T = 5 



(b) T = 10 
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(c) T = 20 



(d) T = 40 



Figure 2: Success probability against minimum gap for a two-qubit system at 
computation times of T = 5, 10, 20 and 40. The Ji have been chosen from the 
Gaussian distribution A/'(0, 1^) for the 500, 000 random problem instances. The 
data points are coloured by | J3I. 



lines correspond to problem Hamiltonians with a degenerate ground state. 
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These lines separate Hf with different ground states: those in the top-right have 
1 11); working clockwise, the others have |10), |00) and |01). 

The second part of Fig. Is] (success probability) is harder to interpret. As 
expected, there is a broad correspondence between lines of low P and those 
of low Amin. The difference, due to second-order corrections, is not presently 
understood. Nor is the interesting observation that Hf with high P are gathered 
around the points of triple degeneracy. 

These plots suggest a projection of a surface onto the (Amin,^) plane; we 
seek to find a suitable parameterisation of the set of Hamiltonians to collapse it 
onto a low-dimensional surface. We find that a plot of P against the minimum 
gap Amin and the position 5* of the gap indeed shows that all points lie close to 
a curved surface (which rises with increasing T). This is understandable, since 





Figure 3: Minimum gap (left) and success probability (right) as a function of 
{Ji^J2) at fixed J3^1.22, T = 5. Each coloured point in the plane is a problem 
Hamiltonian. Observe that points with high success probability (yellow) are 
nested around the points of triple degeneracy. 



those two parameters largely determine the shape of the lowest two energy levels. 
Figure |4] shows a projection of this surface onto the (Amin,^*) plane. Visual 
inspection shows that the colour is to a good approximation a function only of 
position in the plane. Note that the position of the points depends only on the 
Hamiltonian parameters, while the probability depends also on the computation 
time. 




Figure 4: Scatter plot of position 5* of minimum gap against minimum gap 
Amin- The Ji have been chosen from the uniform distribution Z//(— 3,3) for 
100, 000 random problem instances. Points are coloured by the success proba- 
bility P at T = 5. 



4 Larger systems 

We have shown that the relationship between the success probability and Amin is 
not just a case of trivial proportionality for simple two-qubit systems. However, 
it is important to determine whether the interesting structure in this relationship 
remains in larger systems. To determine whether these densely-populated bands 
represent groups of problem instances that have followed similar evolution paths 
for the state vector (e.g. the system remaining mostly in the ground state, then 
being excited at a single avoided crossing), we calculated the average overlap 



with the ground state: 



ds\{0;s\4>{s))\ 



(10) 



The points in Fig. [5] are coloured with respect to this average overlap value, 5^ 
and we can see a smooth graduation across the figures, with the average overlap 
with the ground state increasing with the success probability. The exception to 
the smooth graduation of 6 is the densely populated band where S ^ 1. This 
band must consist of cases with a degenerate or near-degenerate ground state at 
5 = 1, as it includes cases which remain close to the instantaneous ground state 
throughout the majority of the evolution but have a low success probability. 
These results also lends credence to the idea that the structure is closely linked 
to the choice of Hamiltonians. 
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(a) n = 2 



(b) n = 3 
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(c) n = 4 



(d) n = 5 



Figure 5: Distributions of success probability against Amin for two-, three-, four- 
and five-qubit systems over a set of 100, 000 random problem instances, with 
T = 40. The colouring of the points denotes ^, average overlap of the state 
vector with the instantaneous ground state. 



We note that these distributions are reminiscent of the 2D projections of the 
higher-dimensional equilibrium surfaces seen in catastrophe theory [6 . In this 
case success probability, Amin and 6 are all internal variables of the system and 
not independent control parameters, so we are looking at a different situation 
to those usually studied in catastrophe theory. Identifying the nature of this 
surface and the dimensions of the phase space that it exists in is an important 
task, as it could have a major impact on adiabatic algorithm design. At this 
point we can conjecture that the constraint originates from an adiabatic invari- 
ant of the Hamiltonian. We find it strange that, to the best of our knowledge, 
there has been no research on adiabatic invariants of adiabatic quantum com- 
puters. A systematic investigation of adiabatic invariants of quantum computers 



— especially adiabatic and approximately adiabatic computers — could yield 
important information about their behaviour and make a major impact on the 
adiabatic algorithm design. 

5 Conclusion 

We have shown that the relationship between the success probability and the 
minimum ground state gap may not be as straightforward as is often assumed. 
There is a rich structure of distinct sharp edges and densely-populated bands 
in the distribution, particularly in smaller systems. A partial explanation has 
been proposed, whereby this is the projection of a higher-dimensional surface; 
identification of the parameters governing this surface will guide understanding 
of the set of problems amenable to adiabatic quantum computing. We do not 
propose a definitive explanation of the origin of this rich structure: this remains 
an open question. However, we do suggest that there is some evidence that 
some of this structure could arise by our choice of Hamiltonians. We speculate 
that a systematic investigation of adiabatic invariants of quantum computers 

— especially adiabatic and approximately adiabatic computers — could yield 
important information about their behaviour and have a major impact on adi- 
abatic algorithm design. 
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